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A procedure for the instability analysis of three-level multi-span rotor systems 
is described. This procedure is based on a distributed mass-elastic representation 
of the rotor system in several eight-coefficient bearings. Each bearing is 
supported from an elastic foundation on damped, elastic pedestals. The foundation 
is represented as a general distributedmass-efastic structure on discrete supports, 
which may have different stiffnesses and damping properties in the horizontal 
and vertical directions. This system model is suited to studies of instability 
threshold conditions for multi-rotor turbomachines on either massive or flexible 
foundations. The instability condition is found by obtaining the eigenvalues 
of the system determinant, which is obtained by the transfer matrix method from 
the three-level system model. The stability determinant is solved for the lowest 
rotational speed at which the system damping becomes zero in the ccugilex eigenvalue, 
and for the whirl frequency corresponding to the natural frequency of the unstable 
mode. An efficient algorithm for achieving this is described. Application 
of this procedure to a rigid rotor in two damped-elastic bearings and flexible 
supports is described. A second example discusses a flexible rotor with four 
damped-elastic bearings. The third case coqvares the stability of a six-bearing 
300 Mv turbine generator unit , using two different bearing types . These applications 
validate the ccoqauter program and various aspects of the analysis. 


INDECUCnCN 

Problems of rotor instability continue to occur in modern turbomachinery as 
stability limits for speed, power, and flexibility are pressed more closely 
by advanced rotating equipment. Established methods for raising the instability 
threshold speed, such as flexible supports, stable bearing types, etc., are 
sometimes unable to accommodate operational requirements imposed by such designs. 
At the same time, demands for simpler and less costly support structures can 
introduce additional vibration problems to the turbomachine-foundation structure, 
which may further influence the instability threshold speed of the rotor system. 

The purpose of the paper is to describe a general-purpose multi -bearing rotordynamics 
computer code for the calculation of instability whirl threshold speeds for 
three-level rotor systems, of the type shown in figure 1. The rotor has distributed 
mass-elastic properties, and may carry massive disks at the end of each shaft 
section. The usual linear eight-coefficient representation of bearing dynamic 
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properties is employed at each support location, and each bearing may be mounted 
upon a massive pedestal which is flexibly supported with damping from the foundation 
structure. The foundation also has distributed mass-elastic properties with 
allowance for concentrated masses at the ends of each foundation section. The 
foundation is mounted on discrete, damped, flexible supports which attach it 
to ground. The pedestal supports and the foundation supports both attach to 
the foundation, but not necessarily at the same locations: see figure 1. In 
this manner the influence of massive pedestals, and of a flexible , tuned foundation 
on the instability threshold of a multi-bearing rotor may be obtained. 
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Cross sectional area of shaft 

Radial bearing clearance 

Pedestal damping coefficient 

Bearing damping coefficient 

Young's modulus 
Bearing reaction forces 

Cross sectional moment of inertia 

Polar mass moment of Inertia 

Transverse mass moment of inertia 

Pedestal stiffness coefficient 

Bearing stiffness coefficients 

Length of shaft section between station n and n+1 
Bending moment to the left of station n 
Bending moment to the right of station n 
Pedestal mass components 
Mass at rotor station n 
Time 

Shear force to the left of station n 
Shear force to the right of station n 
Bearing static load 

Components of rotor whirl displacement 
Axial coordinate along rotor 
Components of deflected rotor slope 
Angular frequency of rotor 
Critical frequency of rotor 
Frequency ratio = w w /co^ 

Determinant of matrix 
Real part of determinant 
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Imaginary part of determinant 
Threshold frequency, threshold speed 
’Whirl frequency (= U3 ^ * f^^) 

PREVIOUS VDRK 

Studies of the instability properties of multi-bearing rotor systems are 
comparatively rare in the open literature. The earliest rotordynamic analysis 
of such systems appears to date back to Borowicz [1] in 1915, while numerical 
procedures for multi-rotor systems were initiated by Prohl [2] in 1945 for critical 
speeds, and by Koenig [3] in 1961 for multi-rotor system unbalance response 
studies with damped flexible bearings. More advanced multi-bearing rotor system 
procedures and codes were developed by Lund and others, for unbalance response 
and instability analysis of discrete rotor models [4], and for analysis of 
distributed mass-elastic rotors [5] and [6]. Critical speeds of multi-rotor 
systems with flexible bearing were also studied by Crook and Grantham [7 3. 
Other multi-bearing rotordynamic codes were developed by Shapiro and Reddi [8], 
Zorzi and Nelson [9] using the finite element method, and by Gunter [10], Most 
of these studies have included the influence of massive flexible pedestals, 
but there do not appear to have been any published studies in which the influence 
of pedestals and a general mass-elastic foundation on flexible supports is included. 
Most studies have concentrated on the influence of multiple spans on critical 
speeds and unbalance response, or on the influence of massive flexible supports 
on such properties. The only previous multi-bearing flexible support study 
which deals wrath instability threshold speed appears to be that by Lund [6], 
The code upon which Lund's study is based includes pedestal effects, but not 
distributed foundation effects. Lund's study is noteworthy for the experimental 
confirmation of the theoretical threshold speed predictions which it includes. 
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SYSTEM MXEL 

The system model used in this analysis is shown in figure 1. This model allows 
details of the rotor, bearings, pedestals, distributed foundation, and foundation 
supports to be included in the system response results. The number of rotor 
stations and the number of bearing and foundation supports, is limited only 
by the computer space available. Rotor stations and foundation stations are 
independent, and may be included as required writhin the system model. 

The rotor is massive and elastic and its internal damping is small enough to 
be neglected. As the rotor is symmetric about its axis of rotation, it is modeled 
using circular cylindrical sections having mass, elasticity, rotatory and polar ' 
inertia properties distributed along their length. Both transverse bending 
and shear effects contribute to the rotor stiffness. The rotor is referred 
to as the first level of the system. For analysis a rotor is divided into uniform 
shaft sections, and large disks with concentrated mass and inertia properties 
are included at rotor stations at the ends of the shaft sections, i.e., at the 
stations designated 1, 2, 3, etc., in figure 1. For instability studies the 
rotor unbalance is omitted for this model. 

The bearings which support the rotor at specified stations are represented by 
the well-known eight stiffness and dancing coefficient procedure, which includes 
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both direct and cross-coupled relative displacement and relative velocity effects, 
between the journal and the pedestal motions. The bearing properties are therefore 
linear for small displacements, in keeping with the rest of the structural 
components. The pedestals which support the bearings are represented by mass 
properties in the x- and y- directions (no inertia properties), and they are 
supported by elastic, damped members in these directions. This model allows 
two-dimensional transverse pedestal dynamics (without support coupling) to be 
included. The pedestals are the second level of the system model. 

The foundation is a continuous structural member which has different mass-elastic 
properties in the x- 'and y- directions. It is modeled in a similar manner to 
the rotor, using prismatic shear-beam members with distributed properties between 
stations. Discrete mass effects such as casing mass, generator stator, gearbox, 
etc., are incorporated at the end stations. The foundation supports are represented 
by concentrated elastic-damping elements at the ends of the beam sections, as 
required, in a similar manner to the rotor bearing-pedestal supports. The foundation 
supports possess different stiffness and damping properties in the x- and y- 
directions. The foundation is the third level of the system model. 


ECRCES MO BOENIS CN BCHGR STATIONS 

The convention and notation used in this analysis are shown in figure 2, for 
the x-z plane. Similar expressions apply for the y-z plane. The forces and 
moments which act on the n-th rotor station are shown in figure 3. They are 
considered in the following manner. 


Inertia forces. The inertia forces which act on the concentrated mass 
of the rotor station at angular frequency a) are given by: 
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Inertia moments . The inertia moments which act on the concentrated translatory 
and polar inertias at the rotor station at angular frequency w are given 
by: 
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Bearing reactions. Each bearing is represented by eight stiffness and 
damping coefficients in the customary manner. The bearing forces are given 
by: 
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Here x and y are the pedestal displacements, and x and y are the pedestal 
p p * * p ■ P 

velocities. The cross-coupling terms in these expressions couple the bearing 

and pedestal displacements and velocity. The bearing coefficients are 

functions of Sommerfeld number, beating type, and Reynolds number if turbulent 

flow is involved. The influence of bearing moments on the rotor is usually 

small and is neglected in this analysis. 


FORGES AND MOMENTS CN SHAFT SECTIONS 

The analysis procedure for the shaft follows the distributed mass -elastic method 
described by Lund and Orcutt [51, with the difference that the harmonic unbalance 
components are zero in the fourth-order equations of motion. This procedure 
is not repeated here, but is similar to that given in the next section for the 
foundation analysis. 


FORCES AND MOMENTS CXI FOUCATICXI SECTIONS 


A section of the foundation between stations r and r+l is shown in figure 4. 
The foundation analysis is performed by developing a transfer matrix for each 
section, in a similar manner to the rotor analysis. The structure is sub-divided 
into uniform shear-beam sections. All speed-dependent terms in the foundation 
equations are zero, and the foundation is allowed to have different stiffness 
properties in the horizontal and vertical directions. The equations of motion 
for a uniform foundation beam section between end stations in the x- and y- 
directions are 
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To solve these expressions set 


x 


f 


- i co t 

x^. e 


(5) 


Yf = y f 


— iojt 


19 



Substituting gives the expressions 
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The form of these equations is identical, and they are uncoupled because of 
the coordinate symmetry. Tfoere I f I the x- and y- equations must be solved 
separately. To demonstrate the solution procedure for the x- direction, write 
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Substitution in equation 6a gives 
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The solution to this expression is 


x f = cosh X 1 z + sinh X 1 z + cos X 2 z + sin X 2 z (9) 

w * lete D 2 , Dj, and are integration constants to be determined from the 

boundary conditions. The eigenvalues of these expressions are given by 

X 1 L f = XL f { /l + {a X) 4 - (o X) 2 } 1 (l0) 

X 2 L f = XL f { / 1 + (a X ) 4 + (oX)V 

The boundary conditions for the foundation section are 

2 ■ 0i x f ■ x fn* » ' # £n- M f * M f„’ v £ = V £ „ 

A. the other end o£ the section z = L fai z f - . £n+1 . « £ - « £n+1 , M £ . M[ n+1 , 

Vj = V£ n+ j.^ H>on substitution the required transfer matrix expressions for 
the foundation become: 

Whe simplest procedure is to eliminate the distributed rotary inertia effects 
Iy x , Ijy and to account for these terms by using lumped inertia where neces- 
sary arthe end of the foundation section. Thus we will neglect Iy x and Iyy. 
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X f.n+1 b 1n X f,n + L f.n b 3 n G x ,n + h 2n, b 4n M f ,n + h 3n b 7n V f, n 
G f,n+1 ~ 3 V h 2n b 5n L f,n + b 2n G x,n + h 1n b 6n M f,n + h 2n b 4n V f.n 
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2 ... 2 , 


b 2n = (X 2 COsh X 1 L f,n + X 1 COS X 2 L f,n J ‘ (X 1 + X 2> 
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The solution for the foundation in the y- direction is similar to the 
above . 
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TJ-REE-LEVEL SYSTEM 


Details of the rotor-pedestal-foundation models are shown in figure 5. 
the notation shown, the equations governing the harmonic motion of this 
in the x-z plane are: 
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where u> is the rotor exciting frequency. Similar expressions apply in the y-z 
direction. These expressions are used to obtain expressions for the rotor forces 
F zC n> etc., the foundation forces F xcn £, etc., and the pedestal displacements 

x_ , etc., in terms of the displacements of the rotor and the foundation, 
cnp 

The transfer matrix procedure described in the previous section is then used 
to obtain the system transfer matrix. Asstsning that the shear force and bending 
moment are zero at both ends of the structure, the boundary conditions may be 
written: 
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Starting from rotor station 1, the above equations are used to solve for the 
rotor equilibrium as described in reference [111. Each unknown is applied 
separately, by setting x c ^ = 1, with the other displacements zero; then apply 

x gl = 1, others zero etc., through y cl £ = 1, others = zero; y #1 £ = 1, others 

zero. This requires a total of sixteen calculations (eight rotor x,y; eight 
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foundations, x,y). A 16 by 16 matrix is formed from the resulting expressions. 
The determinant of the matrix is the system stability determinant. The lowest 
speed for which this matrix becomes zero is the instability threshold speed. 

In general, there is no particular relationship between the number of stations 
of the rotor and the number of stations of the foundation, or from which locations 
the rotor is supported from pedestals and foundation in relation to the foundation 
support locations themselves. Ihis independent rotor /foundation support situation 
is incorporated within the computer code by allowing a three-level support at 
all rotor and foundation support locations. 


INSTABILITY THRESHOLD ALCERIUM 

The stability determinant contains two unknowns , the rotational speed and 
the whirl frequency ratio fr^. In the general case, the eight dynamic fluid-film 
coefficients are functions of both i) and ft., so that a closed form solution 
is not tractable. The solution is conveniently obtained by iteration as follows. 

Step 1. Obtain points for Ac vs. fr^ and As vs. fr^. 

For a fixed value of oj , the real part Ac and the imaginary part A s of the 
stability determinant are functions of fr ^ and so may be calculated. Several 
points for Ac vs. fr are shown in figure 7(a), and data for A s vs. fr. are 
shown in figure 7(b). 

Step 2. Find zero points of Ac and As by quadratic interpolation. 

As shown in figure 7(a), the program determines two frequency ratios fr^ and 
fr^ where the sign of Ac changes. Then a frequency ratio fr^ which equals 
to (ft. + fr^)/2 is applied for an additional calculation and obtain A c of 
ftj. The three points with f r j, fr^ and f r^ forma quadratic curve which intersects 
the abscissa at f where Ac is zero. Similarly, in figure 7(b), the zero 
point of As is obtained at frequency ratio f 

IT S W 

Step 3. Obtain points for f c (Ac = 0) vs.oj and f s (As = 0) vs.oj. 

Repeat Step 1 and Step 2 with different w until the specified speed range has 
been calculated. Several points for fc(Ac - 0) vs.oj may be obtained, and similarly 
for fs(As =0) vs.oj, as shown in figure 7(c). 

Step 4. Determine threshold speed oj th and corresponding frequency ratio 


Quadratic interpolation is again used. The program determines two speeds 
and w 2 where the sign of (fc - fs) changes. Then a speed oj ? which equals to 
(oj^ -i- oj^)/2 is applied for an additional calculation and obtain fc and fs of 
“3* The three points of fc vs. oj form a quadratic curve, and similarly for 
fs vs. oj . The computer program then determines the intersection point of these 
two curves, i.e.,oj th and frth. oj th is the threshold speed where both Ac and 
As are equal to zero, and frth is the corresponding frequency ratio ofo)^. 

Theoretically, the real parts and imaginary parts of the determinants become 
higher order polynomials in oj and fr, and may have many zero points for each 
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value of (i) . But in practical use, only the zero value encountered at the lowest 
speed is of interest. The computer program organizes the procedure and the 
lowest threshold speed is obtained. 

Other features of this procedure may also be described briefly as follows. 
There are two ways to obtain the real and imaginary parts of the determinant 
of a matrix. Ctae of them is to use a complex number for each determinant 
calculation. But for many digital cooqouters , this will reduce accuracy because 
both real and imaginary parts of the complex number have only about a half of 
the significant digits that the real constant has in the double precision. 
Mother way is to use a real constant for the determinant calculation to obtain 
teal and imaginary parts, and both parts have the same significant digits in 
the double precision. The latter way was used in the computer program developed 
here. Greater accuracy is thereby obtained. 

Because a higher exponent (either positive or negative) may occur on the real 
and imaginary parts of the stability determinant, measures have been taken to 
prevent from computer overflow or underflow in the determinant calculations. 
Furthermore, the eight bearing coefficients are often sensitive to the results 
of the threshold speed, so that they should be checked carefully beforehand. 
A parametric study of variation of bearing coefficients on the result of threshold 
speed has also been performed which confirms the above comments. 


EXMffLE 1 . RIGID ROTOR CN TOD BEARINGS 

The cylindrical rotor shown in figure 8 has been studied by Lund [4] , and is 
supported in damped flexible bearings and pedestals. In this instance this 
system has been mounted upon a continuous foundation mounted on damped flexible 
columns. Details are given in table 1. The threshold speed of the rotor was 
calculated as (a) one-level system (rotor in damped flexible bearings, rigid 
pedestals, and rigid foundation), (b) a two-level system (rotor in the same 
damped flexible bearings, very stiff pedestals , and rigid foundation), (c) as 
a three-level system (rotor, bearings, very stiff pedestals, and very stiff 
foundation and foundation supports), (d) as a three-level system (rotor, bearings, 
very stiff pedestals, and soft foundation and foundation supports) and (e) as 
a two-level system (rotor, bearings and soft pedestals). Computed results for 
these five cases are shown in table 2. The one-level result was checked against 
the one-level result given by Lund. The two-level result (with very stiff pedestals) 
was also checked against the same result . The analysis and results were verified 
through two levels against known data. The effect of very stiff pedestals in 
case (b) did not differ from the case (a) result . The influence of a very stiff 
foundation in case (c) was also found to be small. The influence of a soft 
foundation in case (d) was found to be large, and the effect of soft pedestals 
in case (e) was also found to be large, showing that the effect of pedestals 
and/or foundation can be very important . These results are in agreement for 
a rigid rotor system with a massive foundation. 
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; : n - KX X B XX KXY B X-Y KYX B YX - — KY-Y .... — BY-Y. — 

O. 13E 06 O. 13E 07 0. 62E 06 0. 12E 06 -. 88E 05 O. 12E 06 O. HE 06 0. 18E 06 
BEARING AT STATION 8 

__ . KXX— -BXX — KXY — — BXY KYX- -BYX- KY-Y-- BYY ~ 

O. 48E 06 0. 77E 06 0. 33E 06 O. 16E 06 O. 11E 05 0. I6E 06 0. 12E 06 0. 92E 05 


f r PEDESTAL DATA 

'*’• CTATTHM 


STATION 

MASS-X 

STIFF-X 

DAMP-X 

MASS-Y 

STIFF-Y 

DAMP-Y 

2 - 

0.-100E^17- — 

Q.-l OOE— 11 

. 0.-1 OOE- 03 

0. I00E^17 — 

-0-1 OOE— 1 1 - 

0.-100E-03 

j 8 

0. 1 OOE— 17 

0. 100E 11 

0. 100E 03 

0. 1 OOE- 17 

0. 100E jll 

0. 100E 03 


FOUNDATION SUPPORT. DATA 

STATION STIFF-X DAMP-X STIFF-Y DAMP-Y 

— 2 O. 1 OOE .18 0. 100E jOI... 0.1 OOE .18 _ O. 100E 01 _ 

8 O. 100E 18 O. 100E 01 O. 100E 18 O. 100E 01 


STATIONS WITH INTERCONNECTION BETWEEN ROTOR AND FOUNDATION, FOR ROTOR AS FOLLOW 
2 8 

STATIONS - WITH INTERCONNECTION BETWEEN -ROTOR AND FOUNDATION, FOR . FOUNDATION -AS FOLLOW . __ 

2 8 
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TABLE L- CONCLUDED. 


h: 

FIRST- SPEED - LAST SPEED - - SPEED INCR. _ 

O. S31000E 04 0. 851000E 04 O. 500000E 02 


i’i 
* 1 

I •— FREQUENCY RATIOS - 

V, O. 500000E 00 O. 495000E 00 0. 490000E 00 0. 480000E 00 


TABLE 2.- RESULTS OF EXAMPLE 1 


Case 

System 

Threshold Speed, rpm 

Difference 

Program 

Lund 143 

(a) 

one- level 

8317 

8350 

-0.4% 

<b> 

two-level, stiff pedestals 

83 17* 


-0.4% 

<e) 

three-level, very stiff foun. 

8204** 


-1.7% 

(d) 

three-level, soft foundation 

72X7*** 


-13.6% 

<e> 

two-level , soft pedestals 

7916**** 


-5.2% 


* The two-level system with very stiff pedestals has negligible difference 

in threshold speed from the one-level system. 

** The three-level system with very stiff foundation has a little less 

threshold speed than that of a one-level system. 

*** The three-level system with soft foundation has a much less threshold 

speed than that of one-level or two-level system. This shows where 
the foundation effect is important. 

**** The two-level system with soft pedestals has a somewhat lower threshold 
speed less than that of the one-level system. This also shows pedestals 
effect cannot be ignored. 


BOOTHS 2. FLEXIBLE HOBCR IN EOS EMOTED FLEXIBLE BEARINGS 
CN A FLEXIBLE ECUOOICN 

A four beating totot system given by Lund £4] as an unbalance response example 
is shown in figure 9. Details of the system are given in table 3, The rotor 
was modeled using a 27-station representation and calculated using the eight 
coefficient bearing data. A suitable foundation model was then developed and 
the threshold speed of the rotor was calculated (a) as a one-level system, (b) 
as a two-level system (with very stiff pedestals), (c) as a three-level system 
(with very stiff foundation), .(d) as a three-level system (with soft foundation), 
and (e) as a two-level system (with soft pedestals). Computed results for these 
five cases are given in table 4. The foundation effect was found to be relatively 
small in case (c) and larger in case (d). Similar results were found for the 
pedestals effect in case (b) and (e) . These results further show that in this 
case the effects of foundation exact a significant influence on the whirl threshold 
speed. 
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ORIGINAL PAGE IS 
OF. POOR QUALITY 

TABLE 3.- DETAILS OF THREE-LEVEL, FOUR-BEARING ROTOR SYSTEM EXAMPLE 


TEST: INSTABILITY, 27 STATIONS, 4 BR G. (3 LEVEL) 


INSTABILITY ANALYSIS LEVEL FRQ— RATIO 

2 3 7 


I UN IT IOMEGA ICOEF NLIST ALIST 

0 1 0 0 0. OOE+OO 



DATA 

FOR ROTOR: STATIONS 

BEARINGS 

4 






DATA 

FOR FOUNDATION: STATIONS 

FOUNDATION 

SUPPORTS 






23 

4 






DATA 

OF EACH ROTOR STATION 






»y 

— STA. 

MASS — ROL.- MOM -IRS. -MOM- 

LGTH 

OD. STF_ - 

— — ID-S - 

-OD-MASS 

—ID, M — -YME7 — DENS -SHME7 

■'* 

1 

76. 0 0. OOE+OO 0. OOE+OO 

3. 75 

5. 72 

0. 00 

0.00 

0. 00 3. 130 0. 283 0. 900 


2 

21.4 0. OOE+OO 0. OOE+OO 

4. 09 

5. 62 

0. 00 

0. 00 

0. 00 3. 130 0. 283 0. 900 


-3 

-30. 0 - 0. OOE+OO -0. OOE+OO - 

5.-31 — 

—5. -62 - 

0. 00 - 

0. 00-- — 

~0.-00-3.43G -0. 283- 0—900 — 

i 

4 

125. 0 0. OOE+OO 0. OOE+OO 

6. 62 

7. 36 

0. 00 

0. 00 

0. 00 3. 130 0. 283 0. 900 

”1 

5 

124. 0 0. OOE+OO 0. OOE+OO 

2. 25 

8. 86 

0. 00 

0. 00 

0. 00 3. 130 0. 283 0. 900 




382.-0-0. -OOE+OO- 0 .OOE+OO- 

5 9i 

—11^22- 

0.-00— 

n aa 

0 00 3 130 0 283 0 900 

f 

7 

876. 0 0. OOE+OO 0. OOE+OO 

9. 78 

15. 14 

0. 00 

0. 00 

0. 00 3. 130 0. 283 0. 900 

i 

8 

777.0 0. OOE+OO 0. OOE+OO 

7. 19 

11. 43 

0. 00 

0. 00 

0. 00 3. 130 0. 283 0. 900 


O 

.570r-G— Q,-00E*00— O,-00E+Q0- 

_ g^57 . 

. 9 . 47 

0 . 00 

0 

0 00 3 130 Q 283 X> .900 

rsl 

i 

10 

453.0 0. OOE+OO 0. OOE+OO 

1 7. 19 

9, 41 

0. 00 

6! 00 

0. 00 3. 130 0. 283 0. 900 


11 

235. 0 0. OOE+OO 0. OOE+OO 

5. 00 

8. 99 

0. 00 

0.00 

0. 00 3. 130 0. 283 0. 900 

i,j 

- -12- 

300, 0-0.-00E+0O 0,-00E+0O- 

6.62 — 

- 10. 57 — 

0.-00- 

0. 00- — 

-0. 00 -3. 130 - 0.-283-0. 900 


13 

313. 0 0. OOE+OO 0. OOE+OO 

4. 87 

7. 92 

0. 00 

0. 00 

0. 00 3. 130 0. 283 0. 900 

1 

14 

82. 0 0. OOE+OO 0. OOE+OO 

2. 13 

6. 10 

0. 00 

0. 00 

0. 00 3. 130 0. 283 0. 900 


15- 

— 32, 0- 0, OOE+OO 0. OOE+OO- 

-3.-19 

5.-40- 

0, 00 - 

0,-00 

-0. 00-3. 130 -0.283 0. 900 


16 

36. 7 0. OOE+OO 0. OOE+OO 

7. 42 

4. 49 

0. 00 

0. 00 

O. 00 3. 130 0. 283 0. 900 

i : 1 

17 

28. 2 0. OOE+OO 0. OOE+OO 

8. 65 

3. 25 

0. 00 

0. 00 

0. 00 3. 130 0. 283 0. 900 


—18 

90. 4 O.-OOE+OO O. -OOE+OO- 

12.38 

- - 1. 59- 

0.-00 - 

... — 0 . 00—- 

--0.-00 3.430 0.283 0.900 . ... 

f 

19 

4. 3 0. OOE+OO 0. OOE+OO 

1. 24 

1. 59 

0. 00 

0. 00 

0. 00 3. 130 0. 283 0. 900 

j 

20 

4. 3 0. OOE+OO 0. OOE+OO 

1. 24 

1. 59 

0, 00 

0. 00 

0. 00 3. 130 0. 283 0. 900 

*L. 

2-1- 

78«-5-O,-G0E+OG-Or-0GE+OG~ 

- 4 75 

.,.-4,-49- 

Q, QQ 

0. 00 

n QQ 3 130 Q_ 283 0 .900 

»\: 

_! 

22 

34. 2 0. OOE+OO 0. OOE+OO 

s! 32 

6. 87 

0. 00 

0. 00 

0. 00 3. 130 O. 283 0. 900 


23 

98. 2 0. OOE+OO 0. OOE+OO 

6. 19 

1 1. 48 r ‘ 

0. 00 

0. 00 

0. 00 3. 130 0. 283 0. 900 


—24 

243.-1-0- OOE+OO -O.-OOE+OO- 

6. 51 - 

- 13 - 44 — 

0.-00 - 

0.-00—— 

-0,00 -3,430-0.-283-0^-900 

!J 

25 

330. 0 0. OOE+OO 0. OOE+OO 

6. 21 

11. 42 

0. 00 

0. 00 

0. 00 3. 130 0. 283 0. 900 


’I 26 284. 3 O. OOE+OO 0. OOE+OO 4.33 7. 50 0.00 0.00 0.00 3.130 0.283 0. 900 

M — — 27 63.-7— O »-OOE* GO— O. -OOE+OO— - 0,-00 - —0. 00 0.-00 — - 0. 00— — -0. 00 3. 130 0. 283 - 0. -900- 


- DATA OF - EACH -FOUNDATION- STATION — 


f.j STA. 

MASS 

RSIXY(l) RSIXY(2) 

LGTH 

AXY 

YME7 

DENS 

SHMXE7 

SHMYE7 


0. 0 

O. 24E+11 0. 24E+11 

0. 45 

0. 55E+06 

0. 59E+00 

0. 26E+01 

0. 21E+01 

0. 21E+01 

i i . ... P 

Q O 

A AOCj. 1 1 A AO CTo. 1 1 

A OA A "7 O CTa. A A 

A ®iQP+AA 

A -9Ajr+At 

a oi Pam 

n 01 c 4 .n 1 

i" -3 

63000. 0 

0. 48E+09 0. 4SE+09 

0. 60 

0. 78E+G5 

O. 59E+G0 

0. 26E+01 

0. 21E+01 

0. 21E+01 

U 4 

0. 0 

0. 48E+09 0. 48E+09 

0. 60 

0. 78E+05 

0. 59E+00 

0. 26E+01 

0. 21E+01 

0. 21E+01 

. 5 

-24000.-0 

-0. 20E+10 0.-20E+10 

2.-04 

-0. 16E+06- 

-0.-59E+00- 

-0.-26E+01- 

— 0.-21E+01 

-0, 21E+01 

! 6 

0. 0 

0. 20E+10 0. 20E+10 

2. 04 

0. 16E+06 

0. 59E+00 

0. 26E+01 

0. 21E+01 

0. 21E+01 

7 

62000. 0 

0. 47E+09 0. 47E+09 

0. 65 

0. 76E+05 

0. 59E+00 

0. 26E+01 

0. 21E+01 

0. 21E+01 

: — 8- 

- — 0, 0-0.-47E+09-0-47E+09 — 

— 0. 65 

0. 76E+05 

0. 59E+00 

0. 26E+01 

0. 21E+01 

0. 21E+01 

!' V; 9 

28000. 0 

0. 11E+10 0. 11E+10 

1. 95 

0. 12E+06 

0. 59E+0G 

0. 26E+01 

0. 21E+01 

0. 21E+01 

- 10 

0. 0 

0. 40E+09 0. 40E+09 

1. 96 

0. 71E+05 

0. 59E+00 

0. 26E+01 

0. 21E+01 

0. 21E+01 

r-- - —11—40000.-0 

0, 47E+09 0.- 47E+09— — 

- -0,52 

-0. 76E+OS 

0. 59E+00 

-0. 26E+01 

0. 21E+01 

0. 21E+01 

1 12 

0. 0 

0. 47E+09 0. 47E+09 

0. 52 

0. 76E+05 

0. 59E+00 

0. 26E+01 

0. 21E+01 

0. 21E+01 

:1 13 

148000. 0 

0. 38E+10 0. 38E+10 

2. 85 

0. 22E+06 

0. 59E+00 

0. 26E+01 

0. 21E+01 

0. 21E+01 

! -'L 14- 

O-.-O— Gt-1 2E+08— O.— 12E+08— 

_2.-85 — 0.- 13E+05— 

-0^59E+00- 

-0.-26E+01- 

— 0. 21E+01 

-0.-21E+01 - 


TEST: INSTABILITY# 27 STATIONS, 4 BRQ. (3. LEVEL) 


'T 


| J ! 


15 

16 
. 17 

18 

19 

-20 

21 

22 

-23 


35000. O 

0 . 0 

0-0 

184000. 0 
38000. 0 

3000-0. 

2000 . 0 

O. 0 

0-0 


O. 19E+09 
0, 54E+09 
O— 54E+09- 
0. 19E+09 
0. 31E+08 
0.-31E+08 
0. 99E+07 
0. 26E+09 
O-OOE+OQ 


O. 19E+09 
O. 54E+09 
0.-54E+-09- 
O. 19E+09 
O. 31E+08 
.0, 31E+08- 
0. 99E+07 
O. 26E+09 
.0^ OOE+OO- 


1. 63 

0. 60 
-0-.6G- 

2. 55 
4. 00 

-1—59- 

1. 45 
O. 68 


O. 48E+05 
0. 83E+05 
-O. 83E+05- 
O. 48E+05 
0. 20E+05 
-O- 20E+05 
O. 1 IE +05 
O. 57E+05 


3. 00—0. OOE+OO 


0. S9E+00 
O. 59E+00 
0.-59E+00 
O. 59E+00 
O. 59E+00 
0. 59E+00 
O. 59E+00 
O. 59E+00 
0. 59E+G0- 


0. 26E+01 
0. 26E+01 
- O. 26E+Q1 — 
O. 26E+01 
O. 26E+01 
-0.26E+0t_ 
0. 26E+G1 
0. 26E+01 
0.-26E+01- 


O. 21E+01 
0. 21E+01 
O, 21E+01 
O. 21E+01 
0. 21E+01 
O.-21E+01 
0. 21E+01 
0. 21E+01 
0,21E+01 


O. 21E+01 
0. 21E+01 
-0. 21E+01- 
0. 21E+01 
0. 21E+01 

- 0-21E+01 
0. 21E+01 
0. 21E+01 

- 0. 21E+01- 
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TABLE 3.- CONCLUDED 


_AT_ EACH_OFlFOLLOWING—RQ-TOR_STAIIQti— THERE. -IS - A - BEARING ... . 

3 16 22 27 

AT EACH OF FOLLOWING FOUNDATION STATIONS THERE IS A FOUNDATION SUPPORT 


BEARING DATA .. — .... 

BEARING AT STATION 3 

KXX BXX KXY BXY KYX BYX KYY BYY 

— 0.-1 5E+0.7- 0.. 34E-fcO7.^QOE+06-O^ 13E+07~0^97E+06-0. 1 3E+G7-. 0-1 5E+06 O. 10E+07L-_„„ _ 

BEARING AT STATION 16 

KXX BXX KXY BXY KYX BYX KYY BYY 

0^ 15E+07_CL-34E.-tOX-^aOE-t06_a^t3Ei07_XL-9.7E+06-0^1 3E.+Q7-0- 1 5E+G6~G-10£+07- — 

BEARING AT STATION 22 * r. 

T ' KXX BXX KXY BXY KYX BYX KYY BYY 

... 0. 15E+07 0. 33E+07 10E+07— G— 18E+07— G-4. 2E+07- 0- 18E+G7 -~-%-l 2E+06 0. 18E+Q7 - - — — ^ 

BEARING AT STATION 27 

KXX BXX KXY BXY KYX BYX KYY BYY 

. O. 1 5E+Q7-0. 33E+OZ-t=^10E+07-0._18E+07~0.*12E+07-0. 18E+07 12E+06- 0,~ 18E+G7- 


PEDESTAL 
- STATION 
3 

DATA - 

MASS-X 
0. 100E-G4 

STIFF— X 
0. 100E+14 

DAMP-X 
0. 100E+00 

MASS-Y 
0. 10GE-04 

ST IFF- Y 
0. 100E+14 

DAMP-Y 
0. 100E+00 

16 

0. 100E-04 

O, 100E+14 

0. lOOE+OO 

0. 10GE-04 

0. 100E+14 

0. lOGE+OO 

22 ^ 

0.-1 00E-04 - 

— 0^- 1 00E+ 14- 

■__ -0.- lOOE+OO — 

O-IOOE— 04- - 

_0, 100E+14 - 

o. lOOE+OO 

27 

0. 100E-04 

0. 100E+14 

0: 1Q0E+00 

0. 10GE-G4 

0. 100E+14 

0. 100E+00 


‘ FOUNDAT I ON— SUPPORT— DAT A^ 


STATION 

STIFF-X 

DAMP-X 

STIFF-Y 

DAMP-Y 


5 

0. 100E+33 

. 0. 100E+01 

0, 100E+33 

0. 100E+01 


9 

n i nnp+m n i nnr+m 

Or~l 00E5**33 

0. 100E+33 

-Gr-1 OOE+0 1 — 


13 

0. 10QE+33 

'• G. 100E+01 

0. 100E+G1 

». 

17 

0. 100E+33 

0. 100E+01 

0. 100E+33 

0. 100E+01 

1 


STATIONS WITH INTERCONNECTION BETWEEN ROTOR AND FOUNDATION, FOR ROTOR AS FOLLOW 

... 3 16 22 27 — _ - , — 

STATIONS WITH INTERCONNECTION BETWEEN ROTOR AND FOUNDATION, FOR FOUNDATION AS FOLLOW 
!‘" 4 14 17 21 I 


i 


I 

TEST: INSTABILITY, 27„ STATIONS, 4. BRG. (3 LEVEL) 1 . 


FIRST SPEED V LAST SPEED •*& SPEED INCR. 
O. 900000E+04 O. 100000E+05 O. 200000E+03 


FREQUENCY RATIOS 

0. 500000E+00 O. 495000E+G0 O. 490000E+00 O. 480000E+00 0. 400000E+00 

O. 350000E+00. 0. 300000E+00 — 
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TABLE 4 RESULTS OF EXAMPLE 2 


Case 

System 

'Threshold 
Speed , rpm 

Frequency 
Ratio, fr 

Whirl 
Speed, rpm 

Percentage 
Differ cnee 

(a) 

one- level 

9910 

0.4990 

4945 

0% 

(b) 

two-level 

9910 

0.4990 

4945 

0% 

<c> 

three-level, stiff foun. 

9873 

0.4996 

4943 

-0.24% 

<d> 

three-level, soft foun. 

9625 

0.4998 

4811 

-2.71% 

(e) 

two-level, soft pedestal 

9732 

0.4995 

4861 

-1.7% 


EXAMPLE 3. 300 MW TLIUJINE-GE^JERATOR UNIT 

The turbine-generator rotor shown in f i gure 10 operates at 3600 rpm and is mounted 
in six fluid-film bearings, ranging in size from 200 mm (7.87 in) to 400 m 
(15.75 in) shaft diameter. Bearing oil inlet viscosity is 16 cp at 160oF. 
The bearings are supported on flexible pedestals, mounted upon a stiff foundation, 
which was effectively a two-level system. The rotor model has 15 stations. 
The threshold speed of the system was calculated (a) using two plain cylindrical 
bearings for the generator (i.e., #4 and #5 bearings), and four MPG-Type bearings 
for the other supports, and (b) using six MFG-iype bearings for all supports. 

The results of the stability threshold speed calculations are shown in table 
5. For case (a), the threshold speed is 2200 rpm. The corresponding frequency 
ratio is 0.3518, which corresponds to a whirl speed of 2200 x 0.3677 - 809 rpm. 
This is close to the first flexible critical speed of the generator, which occurs 
at 837 rpm. For case (b), no threshold speed was found below 4000 rpm. These 
results agree with field observations, and demonstrate the capability of MFG-Type 
bearings for preventing rotor whirl instability. 


TABLE 5 RESULTS BOR EXAMPLE 3 


Case 

System 

Threshold 
Speed, rpm 

Frequency 

Ratio 

VShirl Speed 
Nw=Nth x fr rpm 

(a) 

IWo Plain Cylindrical Brgs. 
for Generator, Four MFG Type 
Brgs . for Others 

2200 

0.3677 

809 

(b) 

Six MFG Type Brgs . 

No 

No 

No 


CENCUUSKNS 

1. A general-purpose rotordynamics analysis procedure and computer code has 
been described for the calculation of instability threshold conditions 
for damped flexible rotor-bearing systems mounted in pedestals, and supported 
on a distributed mass-elastic flexible foundation. 
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2. The procedure described has been verified using several established results 
from the open literature, and close correlation between pedestal threshold 
speed and whirl frequencies has been found. 

3. Application of this procedure to a six-bearing, five-rotor, turbine-generator 
system has been described. Results were consistent with field observations 
for this case. 

4. The procedure is suitable for evaluating the influence of pedestals and 
of distributed support foundations on the stability of multi-span rotor 
systems. The transfer matrix procedure is efficient, and double-precision 
accuracy is ensured using the complex eigenvalue algorithm described. 


xxn ommmmr 

The authors are grateful to the Brown Boveri Company of Baden, Switzerland, 
for permission to publish the 300 mw turbine-generator results used in Example 
3. Particular thanks are due to Dr. Anton Roeder , Manager and to Dr. Raimund 
Wohlrab, Senior Engineer, Steam TUrbine Department, for their valuable comments 
and encouragement during this study. 
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Figure 1 Rotor-Bearing-Pedestal Foundation System 
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STATION n STATION (n + l) 


a. End Station and Shaft Nomenclature 



b. Shaft Section Between Rotor Stations n and n+l 


Figure 2 Convention and Nomenclature for End 
Stations and Shaft Sections 


32 





Figure 4 Notation and Convention for Foundation 
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Rotor 



Foundation 


Figure 5 


Three-Level Support System 
for Solution Algorithm 


Figure 6 


Bearing, Pedestal and Foundation 
Model for System 



A c (1) A c change sign between and f r2 




Figure 7 Threshold Speed Determination from the 
Computer 
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Barings 


Pedestals 

PedosUf 
Support* 

Foundation 

Foundation 
Support* 

Figure 8 Rigid Rotor System. Example 1. After Lund [4] 



Figure 9 Flexible Rotor in Four Bearings. After Lund [4] 




,p L “ Generator Exciter 

Figure 10 300 Mw Turbine-Generator Rotor 

System. Five Rotors; Six Bearings 





